
f1=new TFile("histRecTHnIncl_LHC10-d-e_TPC0TOF0_ETA0_JPTall_0-100/histRecTHnIncl_LHC10-d-e_TPC0TOF0_ETA0_JPTall_0-100.root")
htpc0=horiginalyield
f2=new TFile("histRecTHnIncl_LHC10-d-e_TPC1TOF0_ETA0_JPTall_0-100/histRecTHnIncl_LHC10-d-e_TPC1TOF0_ETA0_JPTall_0-100.root")
htpc1=(TH1D*)f2->Get("horiginalyield")
htpc1->Divide(htpc1,htpc0,1,1,"B")
gStyle->SetOptStat("")
style::ToNaturalScale(htpc1)
htpc1->Draw()
htpc1->SetXTitle("p (GeV/c)")
htpc1->SetYTitle("raw yield ratio (Geo./Basic)")
htpc1->SetTitle("LHC10 d+e, |#eta|<0.2, inclusive")
gPad->SetLogx(1)
c1->Print("GeoEff_inclusive_ETA0.eps")
c1->Print("GeoEff_inclusive_ETA0.C")

f1=new TFile("histRecTHnIncl_LHC10-d-e_TPC0TOF0_ETA3_JPTall_0-100/histRecTHnIncl_LHC10-d-e_TPC0TOF0_ETA3_JPTall_0-100.root")
htpc0=horiginalyield
f2=new TFile("histRecTHnIncl_LHC10-d-e_TPC1TOF0_ETA3_JPTall_0-100/histRecTHnIncl_LHC10-d-e_TPC1TOF0_ETA3_JPTall_0-100.root")
htpc1=(TH1D*)f2->Get("horiginalyield")
htpc1->Divide(htpc1,htpc0,1,1,"B")
gStyle->SetOptStat("")
style::ToNaturalScale(htpc1)
htpc1->Draw()
htpc1->SetXTitle("p (GeV/c)")
htpc1->SetYTitle("raw yield ratio (Geo./Basic)")
htpc1->SetTitle("LHC10 d+e, 0.6<|#eta|<0.8, inclusive")
gPad->SetLogx(1)
c1->Print("GeoEff_inclusive_ETA3.eps")
c1->Print("GeoEff_inclusive_ETA3.C")
